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Neutron scattering is used to study NiCl2-2ccBr2x-4SC(NH2)2, x = 0.06, a bond-disordered modi¬ 
fication of the well-known gapped S = 1 antiferromagnetic quantum spin system NiCl2-4SC(NH2)2. 

The magnetic excitation spectrum throughout the Brillouin zone is mapped out at T = 60 mK 
using high-resolution time-of-flight spectroscopy. It is found that the dispersion of spin excitation 
is renormalized, as compared to that in the parent compound. The lifetime of excitations near the 
bottom of the band is substantially decreased. No localized states are found below the gap energy 
A ~ 0.2 meV. At the same time, localized zero wave vector states are detected above the top of 
the band. The results are consistent with a more or less continuous random distribution of bond 
strengths, and a discrete, possibly bimodal, distribution of single-ion anisotropies in the disordered 
material. 

PACS numbers: 75.10.Kt,75.10.Jm,75.40.Gb,78.70.Nx 


I. INTRODUCTION 


Quantum magnets with disorder and randomness are 
presently attracting a great deal of attention.EJ® They 
often exhibit behavior that is qualitatively different from 
that of their disorder-free counterparts. There are ac¬ 
tually several ways to introduce randomness in a mag¬ 
netic material. The effect of site disorder (or site di¬ 
lution) is rather well understood. Randomly removing 
spins has particularly severe consequences for gapped 
quantum paramagnets. Their non-magnetic ground state 
is destroyed, and magnetic order often emerges at fi¬ 
nite temperatures through the formation of correlated 
clusters around the impuritics.l^^® In contrast, the ef¬ 
fect of random exchange interactions, or bond disorder, 
is more subtle. In systems with a spin gap the ground 
state may survive. However, materials close to a mag¬ 
netic quantum-critical point (QCP) are more drastically 
affected. Indeed, near QCP even a small variation in 
the Hamiltonian parameters niay lead to the qualitative 
change of the ground state.^^®^ Recently, VojtsP^ has 
numerically investigated the effect of bond disorder on 
the dynamic properties of a gapped quantum magnet 
in the vicinity of the z = 1 QCP.^ Among the pre¬ 
dictions are some intriguing disorder-induced features, 
such as localized states inside the gap, or the possi¬ 
bility of “weak ordering” Unfortunately, despite the 
rapid increase in the number of organometallic quan¬ 
tum magnets,^ and the ease of achieving bond disor¬ 
der via chemical substitution on nonmagnetic sites in¬ 
volved in super exchange there is presently a short¬ 
age of suitable model compounds to test the predic¬ 
tions of Ref. n Most of the known gapped quantum 
magnets (e^. TlCuCls!^ IPA-CuClaP PHCCji^ 
Sul-Cu 2 Cl 4 PCu(Qnx)ClP>^ exhibit a “blue shift” of 
magnons in the presence of bond disorder. Thus, chemi¬ 
cal modification increases the spin gap and pushes these 


systems away from the z = 1 QCP. 

One known exception is the extensively studied S = 1 
gapped quantum magnet dichlorotetrakis-thiourea nickel 
NiCl2'4SC(NH2)2 (commonly abbreviated as DTN).1^^H^ 
In contr ast to integer-spin Heisenberg magnets such as 
CsNiClaP^ the gap in DTN is due to huge planar 
single-ion anisotropy. The three-dimensional interac- 
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FIG. 1. (Color online) A sketch of the 

NiCl2-2ccBr2cc-4SC(NH2)2 structure. The surrounding 
thiourea units are shown only for one of the Ni^^ ions. The 
main magnetic interactions (see Hamiltonian 0 ) for S = 1 
nickel ions are labeled. Two non-equivalent chlorine positions 
are marked as 1 and 2. The scattering plane for the present 
neutron experiment is shown on the right, in both real (top) 
and reciprocal (bottom) space representations. 
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tions in DTN are also stro ng en ough to make the Hal¬ 
dane gap physics irrelevant Bond disorder is intro¬ 
duced into DTN by a random Br substitution on the 
Cl site. This act ually produces a “red shift” (decrease) 
of the spin gapl ^^ l ^^ * In this respect, Br-substituted 
DTN, NiCl2-2a:Br2a;'4SC(NH2)2 (below abbreviated as 
DTNX), appears to be the most promising candidate for 
an experimental realization of the physics discussed in 

Ref.im 

The present work is a study of spin dynamics in DTNX 
with X = 6% Br substitution. High-resolution inelastic 
neutron scattering experiments provide a detailed pic¬ 
ture of magnetic excitations throughout the Brillouin 
zone. The measured spectra exhibit certain key differ¬ 
ences compared to the parent material: an increase of 
the bandwidth, a reduction of excitation lifetimes, and 
new localized states at high energies. We discuss the re¬ 
lation between these experimental results and the numer¬ 
ical predictions of Ref. [121 In addition, we touch on the 
relevance of our findings to the previous th ermod ynamic 
studies of DTNX in applied magnetic fields 

The paper is organized as follows: In Sec. [H] we re¬ 
view the structure and magnetic properties of our target 
material, and describe the experimental setup and tech¬ 
niques; a brief overview of the collected data is given in 
Sec. |HI[ followed by detailed analysis and discussion in 
Sec. |IV| Sec. [V| summarizes our conclusions and outlines 
the future research directions; some more technical issues 
are laid out in the appendixes. 


II. MATERIAL AND EXPERIMENTAL 
DETAILS 

A. The parent and disordered compounds 

Disorder-free DTN NiCl2-4SC(NH2)2 is an 
organometallic magnet belonging to a highly sym¬ 
metric 14 tetragonal space group .1^ Magnetic S = 1 
Ni2+ ions occupy the positions in the body-centered 
tetragonal lattice (see Fig. with a = 9.558 A and 
c = 8.981 A. Due to this body-centered arrangement, 
the magnetic system can be seen as a combination of 
two interpenetrating tetragonal lattices that are totally 
decoupled at the mean-field level. Within each of these 
subsystems the magnetic properties are captured by the 
following Hamiltonian: 




(1) 


where (i,jf) is the summation over the nearest neigh¬ 
bors in the tetragonal (a, a) plane, and index n refers to 
the summation along the c axis. The magnetic properties 
are dominated by the strong single-ion anisotropy of easy- 
plane type D 0.78 meV. The next term in the magnetic 
hierarchy is the nearest-neighbor Heisenberg exchange 


interaction Jc — 0.141 meV within the chains running 
along the high-symmetry axis. The least important inter¬ 
action is the interchain coupling Ja — 0.014 meV, which 
is an order of magnitude weaker than Jc-^ As a result 
of the huge anisotropy, the ground state of DTN is a 
quantum-disordered XY paramagnet with a spin-singlet 
nonmagnetic ground state and a spin gap A 0.3 meV. 

It has been previously demonstrated that Br substi¬ 
tution has a very slight effect on the lattice parameters 
and does not affect the symmetry.l^ The bromine ions 
are site selective: They occupy Cl-1 positions in the lat¬ 
tice as shown in Fig. Chemical substitution decreases 
the critical field Hd of mag netic ordering in DTNX, in¬ 
dicating a reduced spin gap j^^ l ^^ * This gap reduction is 
significant: Already at x = 13% bromine concentration 
A 0.1 meV is almost three times smaller than in the 
parent compound. 


B. Quantities measured 

Magnetic inelastic neutron scattering almost directly 
probes the dynamic spin structure factor. 

5“^(Q,w) = j (5“(o,0)S'^(r,i))^^. 

( 2 ) 

The actual neutron intensity measured in experiments 
is proportional to the differential cross section, and for 
unpolarized neutrons is given by® 


a ^ ^ / 

(3) 

Here Q and huj the momentum and energy transfers, 
respectively. F‘^{Q) is the magnetic form factor for Ni^+, 
known from numerical calculations.® No absolute nor¬ 
malization of the data was performed in this study. 


C. Experimental setup 

In the present work we used large (mass r\j 1 g) fully 
deuterated single crystals of DTNX with 6% Br substi¬ 
tution. They were grown fr om the so lution by the same 
method as in previous worksThe measurements 
were carried out on the high resolution cold neutron 
time-of-fiight (TOF) spectrometer IN5 at Institut Laue- 
Langevin.® Sample environment was a ^He-^He dilution 
cryostat. All data were collected on two co-aligned sin¬ 
gle crystals at T = 60 mK, which is much lower than all 
the relevant energy scales of DTN. The principal exper¬ 
imental scattering plane was (1, —1,0), providing access 
to scattering vectors of type (h, h,/) (see Fig. [^. Two 
data sets were collected, using incident neutron energies 
E[ = 2.26 and 1.17 meV, respectively. A Gaussian fit to 













3 


the elastic incoherent sc atteri ng provided an estimate of 
the energy resolution: Vs\n2a ^ 38 jneV and 16 jneV full 
width at half height, respectively. In the high-resolution 
setup, incoherent elastic scattering has virtually no effect 
on the data collected above hw = 25 /ieV energy trans¬ 
fer. For each incident energy, the sample was rotated 
step-wise through the data collection, to fully cover the 
first Brillouin zone. 


III. DATA OVERVIEW 

The bulk of the neutron TOF data collected in 6% 
DTNX at T = 60 mK is visualized in Fig. These are 
a series of false color plots of the measured inelastic in¬ 
tensity, presented as two-dimensional momentum-energy 
cuts along high-symmetry directions. The data below 
and above 0.7 meV energy transfer were taken with neu¬ 
tron incident energies Ei = 1.17 meV and Ei = 2.26 meV, 
respectively. We stress that the data are plotted “as is”, 
without any background subtraction.!^ 

The measured spectrum of 6% DTNX is qualitatively 
similar to that of disorder-free DTN.I^It is dominated by 
a well-defined magnon mode with a dispersion minimum 
at the antiferromagnetic zone center O = (1/2,1/2,1/2). 
Two saddle points are located at M = (0, 0,1/2) and A = 
(1/2,1/2,0), correspondingly. The highest energy of the 
single-magnon excitation is found at X = (0,0,0). The 
main difference with the parent material is a much re¬ 
duced spin gap A 0.2 meV, compared to A 0.3 meV 
in DTN. 

At lower energies, the excitations in DTNX appear less 
sharp than at high energy transfers. Notably, there are 
no additional excited states below the gap energy at the 
O point or any other region of the Brillouin zone down to 
25 /ieV. Finally, we note a subtle drop of intensity around 
0.4 meV. The latter is found in data sets with both neu¬ 
tron incident energies, and hence is unlikely to be of 
instrumental origin. Although it resembles an avoided 
crossing in the false-color plot shown, it is actually al¬ 
ways observed at the same energy transfer, irrespective 
of wave vector. Its origin remains unclear. 


IV. ANALYSIS AND DISCUSSION 
A. Dispersive excitations 

Since we lack a suitable microscopic model to glob¬ 
ally describe the spectrum in a quantum magnet with 
disorder, we chose a more empirical approach to data 
analysis. We quantitatively analyze the E^ = 2.26 meV 
dataset, covering the full magnetic excitation band. The 
measured data were broken up into a series of individ¬ 
ual constant-Q cuts. In each such cut, the scattering 
is a well-defined peak that we approximated by a Voigt 
profile: 


/(Q, u) = SqF^Q) (^2 - Vihco -hwQ,a, Tq). 

(4) 

The Gaussian width a represents experimental energy 
resolution (the values quo ted a bove, with the energy 
transfer correction included^^^^ . For each cut, the pa¬ 
rameters of this model are as follows. Hujq is the position 
of the peak that we associate with the single-magnon 
energy. The Lorentian width Fq of the Voigt function 
represents the intrinsic magnon line width and, poten¬ 
tially, wave vector resolution (“focusing”) effects. Sq 
is the peak’s integrated intensity. Due to the planar 
nature of the ground state, we assumed that the ob¬ 
served magnon corresponds to transverse spin fluctua¬ 
tions 5^^(Q,co’) = 5^^(Q,u;). The polarization factors 
in the expression above are chosen accordingly, with Q± 
being the momentum transfer in the (a, a) plane. Exam¬ 
ples of fits to individual cuts can be found in Fig. 

The thus obtained magnon dispersion relation is plot¬ 
ted in symbols in the upper panel of Fig. These data 
were further analyzed within the random-phase approx¬ 
imation (RPA),!^ treating the Heisenberg exchange as a 
perturbation to decoupled S = 1 single ions. The RPA 
dispersion relation is: 

hwQ = ^D^+4DjiQ). (5) 

Here 7(Q) = ^ Jr cos(Qr) is essentially the Fourier 

r 

transform Heisenberg exchange interactions in the sys¬ 
tem. For our fit, we consider not only the parame¬ 
ters D, Jc and Ja, but also two possible perturbations 
to the Hamiltonian 0 : diagonal exchange between the 
tetragonal sublattices Jd and the next-nearest neighbor 
exchange along the c direction Jc 2 • The reasons to intro¬ 
duce these perturbations are discussed in Appendix 
Fitting the RPA dispersion relation to the data in FigT^ 
(top panel) yields the parameter values summarized in 
Table[T[ The fit itseif is represented by the soiid iine. For a 
direct comparison, we aiso quote Hamiitonian parameters 
for disorder-free DTN from the previous studie^^ and 
piot the corresponding dispersion reiation in a dashed 
iine. 

As the anaiysis shows, the gap softening in DTNX is 
primariiy due to the increase of the bandwidth. Even as 
the anisotropy becomes stronger with the Br substitu¬ 
tion, the ratio D/Jc actuaily decreases from 5.5 in pure 
DTN to 5.1 — 5.4 in 6% DTNX, depending on the terms 
inciuded in the modei Hamiitonian. The ratio Ja/Jc 
is stiii around 0.1 within the uncertainty iimit. The 
RPA dispersion 0 can be modified to account for quan¬ 
tum renormaiization of the bare Hamiitonian parame¬ 
ters. This correction, described in Appendix [B| ch anges 
the D/Jc ratio to 4.7 in the parent materiai.!^^^ In 6% 
DTNX, as found from the present experiment, it wiii be 
between 4.06 and 4.36. The ratio JajJc — OT is un- 
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FIG. 2. (Color online) False color plot of the neutron scattering intensity measured in 6% DTNX at T = 60 mK. The data 
are presented as a multi-panel cut along the trajectory shown in the inset. The data collected with Ei = 2.26 and 1.17 meV 
incident energy are combined in each panel. The integration range in Q transverse to the cut directions is ±0.02 reciprocal 
lattice units. 


DTN 



6% DTNX 



Zapf 

Present 

study 

Present study Present 

study 

Present study 

ei al. (2006j^ 

Eq. n 

only 

Eq. n±JcZ Eq. H 

•±«/c2 

Eq. (fl|)±e/ci Jc 2 


D 

0.780(3) meV 

0.792(3) meV 

0.812(2) meV 

0.805(2) meV 

0.807(2) meV 

Jc 

0.141(3) meV 

0.155(1) meV 

0.158(1) meV 

0.149(1) meV 

0.150(1) meV 

Ja 

0.014(1) meV 

0.0158(3) meV 

0.0161(2) meV 

0.0155(2) meV 

0.0157(2) meV 

Jd 

- 

- 

0.0059(5) meV 

- 

0.0060(4) meV 

Jc2 

- 

- 

- 

-0.0095(7) meV 

-0.0096(7) meV 


D/Jc 5.5(1) 

5.11(3) 

5.14(3) 

5.40(3) 

5.38(3) 

Ja/Jc 0.098(7) 

0.102(2) 

0.102(2) 

0.104(2) 

0.105(2) 


TABLE I. Results of fitting the RPA dispersion relation ® to the measured excitation energies in 6% DTNX. Additional 
perturbations to the basic Hamiltonian are considered (see Appendix [A| for more details). The results of an analogous 
analysis for disorder-free DTN from Ref. |27| are also shown. 


changed. Thus the reduced D/Jc ratio is the key ingre¬ 
dient of the excitations “red shift” in DTNX. 

The intensities Sq obtained in fits to individual cuts 
are plotted in symbols in the lower panel of Fig.|^ In the 
RPA, Sq is simply proportional to I/JEjq. As shown by 
the solid line, this expectation is in a reasonable agree¬ 
ment with the data. The remaining discrepancies might 
be attributed to the possible contribution of longitudinal 
fluctuations, not included in our approach. 


B. Approaching criticality 

Due to a random distribution of Br substitutes in the 
sample, the parameters of the microscopic Hamiltonian, 
including exchange constants and anisotropy, will them¬ 
selves vary from one unit cell to the next in a random 
manner. This said, it stands to reason that the values 


obtained from analyzing magnon dispersion curves cor¬ 
respond to the average values of these parameters. The 
main mechanism leading to a decrease of the spin gap 
and driving DTNX closer to the QCP is a then a steadily 
decreasing D/Jc ratio. 

Fortunately, the pha se diagr am of Hamiltonian is 
well known numerically, F^F9 EQ] ^^d shown in Fig. |5 For 
small anisotropy and almost isolated chains it includes 
a gapped topological Haldane phase. Sufficiently strong 
interchain interactions restore XY-like long range mag¬ 
netic order. However, the system is again in a gapped 
non-magnetic “single ion” state for large anisotropy. The 
values of Hamiltonian parameters determined for the par¬ 
ent DTN compound,^ allow us to place it in the latter 
region of the phase diagram (see Fig. |^. As mentioned 
above, Br substitution does not change the Ja/Jc ratio 
significantly. Instead, by decreasing DjJ^^ the system 
is driven left on the phase diagram, towards the line of 
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FIG. 3. (Color online) Constant-Q cuts at the center (O) and 
at the boundary (A) of antiferromagnetic Brillouin zone of 6% 
DTNX. Squares and circles correspond lo Ei — 1.17 meV and 
Ei = 2.26 meV, respectively. Voigt hts @ are shown by solid 
lines. The shaded areas represent the experimental energy 
resolution. 


long-range XY ordering. The inset of Fig.j^show the en¬ 
ergy gap A in DTNX as a function of Br co ntent, as de¬ 
duced from the current and previous studies.^^^^^^^ The 
decrease is roughly linear, and the gap may be expected 
to close around 20% Br content. At present it is not 




FIG. 4. (Color online) Top panel: Magnon dispersion relation 
measured in 6% DTNX, as deduced from Voigt hts to individ¬ 
ual constant-Q cuts of = 2.26 meV dataset. The solid line 
is a ht to the RPA prediction of Eq. with the parameters 
summarized in the last column of Table [J The dashed line 
is the reference magnon dispersion in disorder-free DTN fol¬ 
lowing Ref. 1271 Bottom panel: Integrated intensity corrected 
for the magnetic form factor and polarization factor, <Sq, as 
determined in hts to the same individual cuts. The solid line 
is the RPA result, which is simply the inverse of Houq. 


clear whether DTNX samples with such high Br concen¬ 
trations are stable. 


C. Magnon lifetimes 

Beyond a simple change of average Hamiltonian pa¬ 
rameters, several key features of the spectrum are to be 
attributed to disorder, z.e., to a microscopic random vari¬ 
ation of these parameters in the sample. As shown in 
Fig-i the excitations at the antiferromagnetic zone cen¬ 
ter O have a signihcant energy width beyond the reso¬ 
lution of the instrument. This is further emphasized by 
comparing data obtained with different neutron incident 
energy. For E^ = 1.17 meV the energy resolution is much 
sharper than for Ei = 2.26 meV, but the obtained broad 
peak at Q = (1/2,1/2,1/2) is similar in both cases, with 
effective Lorentian linewidth F = 34±4 /ieV. Although at 
high energies the excitations are sharper, their linewidth 
is still pronounced. Fitted linewidth at, for instance, 
point A is F = 22 ± 3 meV. As the resolution function of 
a time-of-flight sp ectrom eter sharpens with the increase 
in energy transferthe linewidth F still exceeds the 
estimated instrument resolution a at Jiw 1 meV almost 
twice. 

Wavevector resolution effects alone can not account 
for the observed increase of linewidth. A conservative 
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FIG. 5. (Color online) Ground states of Hamiltonian 0 - 
based upon the results of Refs. [49] and [SO] The arrow indicates 
the approximate trajectory of DTNX upon Br substitution. 
The inset shows the energy gap in DTNX as a function of 
bromine content. D ata poi nts are the result of the present 
and previous studies.^ ^ * ^^ * ^^ * The dashed line is a guide for the 
eye. 



FIG. 6. (Color online) False color map of measured intrinsic 
Lorentian magnon linewidth Fq as found from Voigt fits to 
individual constant-Q cuts (Ei = 2.26 meV dataset). Regions 
of reciprocal space with no data are shown hatched. 


estimate of Q-resolution contribution to the peak width 
in energy cuts can be done as ctq < \JiWQ-^sQ — In 



(0, 0, l) (r.I.u.) 


FIG. 7. (Color online) High contrast false color map of 
/(Q,ct;) for Q along c* (here Ei = 2.26 meV). The integration 
range in transverse direction is the entire Brillouin zone. Note 
the dispersionless (local) excitations just above the top of the 
magnon band denoted by arrows. 


our case, for the antiferromagnetic zone center this gives 
only ctq 5 /ieV, which is much less than the observed 
broadening. 

The contour map of measured intrinsic line width Fq is 
shown in Fig. One can see that Fq varies significantly 
along the tetragonal axis, but has a less pronounced vari¬ 
ation in a transverse direction. The broadest excita¬ 
tions are found around the antiferromagnetic zone center 
where the dispersion is a minimum. This is in a strong 
contrast with the picture of excitation broadening previ¬ 
ously observed in another bond-disordered gapped anti- 
ferromagnet PHCX.I^ In the latter case, magnon damp¬ 
ing is due to scattering on isolated (discrete) impurities. 
As a result, it roughly scales with the single-magnon den¬ 
sity of states, and is a maximum at the band top. In 
DTNX we find exactly the opposite. As will be argued 
below, the observed behavior of Fq is consistent with a 
smooth continuous distribution of Heisenberg exchange 
strengths in the material. 


D. Localized states 

Figure shows a high-contrast false color plot of 
/(Q, cj) integrated over the whole Brillouin zone and pro¬ 
jected onto a c* direction. Note that here we use a lin¬ 
ear intensity scale, compared to the logarithmic scale 
used above in Fig. In Fig. the magnon band 
looks completely “overexposed”. There are clearly no 
additional states visible inside the gap. However, addi- 
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Both D and J are increased 
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localized state 
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FIG. 8. (Color online) A cartoon highlighting the main effects 
of various kinds of disorder on the spin excitation spectrum. 
The solid line is the magnon dispersion in the absence of dis¬ 
order. The blurred curve is the spectrum of the disordered 
system. 


tional states are found just above the top of the band at 
liw 1.25 meV. These are local excitations, as they show 
no dispersion. 

Given the low intensity of the feature in 6% DTNX, 
how reliable is this observation? An important argument 
here is that the feature at 1.25 meV has an intensity dis¬ 
tribution, perfectly matching the sample’s reciprocal lat¬ 
tice. The intensity is cosine modulated along [0 0 1] and 
constant along [1 10]. This is nontypical for the spurious 
features usually having no Q-structure.^^ As the numer¬ 
ical calculations show, localized feature of this kin d can 
actually be expected in a bond-disordered systenP^^^ 
(this will be discussed in more details below). Indi¬ 
rect evidence of high-energy states in DTNX also comes 
from the experimentally observed pre-saturation “pseu¬ 
doplateau” in magnet iz at ion.The width of the plateau 
AH 1.5 T and the observed energy separation between 
the band top and the localized state hAoo ^0.2 meV are 
in a rough agreement. As these states are located at the 
regions of reciprocal space, equivalent to Q = 0 momen¬ 
tum transfer, they should also be observable by such tech¬ 
niques as electron spin resonance and THz spectroscopy. 
A very recent THz spectroscopy experiment does show 
the presence of an additional feature at approximately 
the same energy in DTNX and its absence in the parent 
material.l^ 


E. Disorder analysis 

To understand the emergence of local excitations and 
other features of the observed spectrum, we can establish 
a crude qualitative “mapping” between the Hamiltonians 
of DTN 0 and that of the dimer modeled studied in 
Ref. [121 The latter has two parameters: the intradimer 


exchange constant J and the interdimer coupling JC. The 
dimer strength J primarily determines the gap and in¬ 
terdimer exchange 1C determines the bandwidth. The 
increase of the former drives the system away, and of 
the latter — closer to QCR In this sense, our S = 1 
individual magnetic sites with single-ion anisotropy can 
be seen as spin-gap objects analogous to the S = 1/2 
dimers. The parameter J is then naturally mapped to 
the single-ion anisotropy D and the critical coupling 1C 
corresponds to the Heisenberg exchange J. Similar map¬ 
ping between dimerized and single-ion anisotropy sys¬ 
tems has also been used in a recent theoretical study 
of Utesov et al^^ 

The numerical study Ref. [12] predicts the presence of 
in-gap states in two cases of discrete disorder distribu¬ 
tion: a small fraction of isolated sites having J < {J) 
(in our correspondence with DTN D < (D))P^or a small 
fraction of sites with 1C > {1C) (in our case J > (J)).l^ 
Our analysis of the magnon spectrum in DTNX shows 
that anisotropy D is effectively increased by Br substi¬ 
tution, so the former scenario is clearly not applicable. 
Thus, not having any in-gap bound states in DTNX may 
indicate that the disorder of the exchange constant J 
is either weak or non-discrete. At the same time, as 
captured by the simulations of Ref. |T2] (Fig. 5 therein), 
a small number of sites with a singular large value of 
J > {J) (i. e. D > (D)) produces a localized state 
near the top of the magnon band. This is totally in line 
with our observation of high-energy local excitations in 
DTNX. The underlying assumption is that the anisotropy 
distribution is discrete^ with just a few Ni^+ ions having 
a substantially increased D-term. The simulations also 
predict a broadening of the main magnon branch, in the 
case of a continuous broad distributions of either 1C and 
J . The former broadens the magnons evenly in the entire 
Brillouin zone (Ref. [12] Supplementary Material, Fig. S4), 
while the latter mostly affects magnons at the bottom of 
the band (Fig. S5). Comparing it to our observations, we 
may again guess that the exchange constants in DTNX 
show a rather broad continuous distribution. 

Certainly, the above discussion is based on a rather 
tenuous analogy between two very different Hamiltoni¬ 
ans. However, if we accept this qualitative correspon¬ 
dence, the following picture emerges. The average values 
of both {D) and (J) are increased with the Br substitu¬ 
tion, modifying the magnon dispersion relation accord¬ 
ingly. The anisotropy is decreased on only a few sites, 
probably in direct proximity of the Br substitutes. The 
exchange constants, however, have a broad statistical dis¬ 
tribution around the average value. Schematically, this 
scenario is borne out in Fig. A broad distribution of 
exchange constants in DTNX is not unexpected. It re¬ 
sults from each Br substitute affecting a large number 
of bonds. Each such size-mismatched Br defect creates a 
strain field in the crystal, which, in turn, affects bond an¬ 
gles and thereby superexchange interactions. The strain 
field falls off as ^ r~^ (Ref. [56|) . and in a soft material 
like DTNX is expected to be rather long-range. 
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FIG. 9. (Color online) Symbols: squared dispersion 
along the (/i,/i,0) and (h,/i, 1) directions {Ei = 2.26 meV). 
The dashed line is the dispersion model for pure DTN used 
by Zapf et a/.^. The solid line is best fit with Jd perturbation 
included, as in the last column of Table 


V. CONCLUSIONS 

In summary, Br substitution has a profound effect on 
the spin dynamics of DTN, even at rather low concentra¬ 
tions. On the one hand, both anisotropy and exchange in¬ 
teractions are, on the average, increased. The ratio D/Jc 
actually decreases, reducing the spin gap and driving the 
system closer to the QCP. Simultaneously, new features 
emerge due to disorder. Magnon lifetimes are shortened, 
predominantly at the bottom of the band. Somewhat 
counter-intuitively, localized states appear not inside the 
spin gap, but just above the top of the magnon band. 
This behavior can be explained by a hand-wavin g a nal- 
ogy with disorder in the Heisenberg-dimer model.^ We 
hope that our experiments will stimulate a study of ran¬ 
domness in the anisotropic single-ion Hamiltonian appro¬ 
priate for DTN. 
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h (r. I. u.) 


FIG. 10. (Color online) Contour map of weighted peak in¬ 
tensity Sq ■ hwQ in the plane {Ei = 2.26 meV). Note 

the half-period modulation along c*, i. e. the peaks at both 
integer and half-integer values of 1. 


the ESR experiment by Zvyagin et and confirmed 
in neutron scattering investigation by Tsirulin et al^ 
The introduction of Jd leads to additional double-period 
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icity can be indeed spotted in the data. Figure shows 
a comparison of dispersion curves (actually, dispersion 
squared) measured along the (/i,/i,0) and (/i,/i, 1) di¬ 
rections (symbols). The difference between the other¬ 
wise equivalent reciprocal-space directions is naturally 
explained by a non-zero Jd^ as shown by the solid line 
which corresponds to the best-fit value of Jd- For com¬ 
parison, the dashed line is the proposed DTN dispersion 
from Ref. [23 where Jd is not included. 
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Appendix A: Additional interactions 

Here we describe some features in our inelastic data 
that point to the presence of additional magnetic inter¬ 
actions, not included in the Hamiltonian Q. The first 
additional interaction is the so-called “diagonal” Heisen¬ 
berg exchange between the two tetragonal sublattices. 
The presence of this exchange Jd was already noted in 


The second additional term that we consider is the 
next-nearest neighbor exchange coupling Jc 2 along the 
chain direction. The evidence for this term comes from 
the very specific pattern found in the intensity distribu¬ 
tion. The product • JEjq shown in Fig. [Tol exhibits 
a half-period modulation along the c* direction. This 
implies the presence of additional interaction with the 
translation vector r = 2c. Additionally, J 2 c manifests it¬ 
self in a slight “blunting” of the cos-like dispersion along 
the c axis, near the top of the band. That, too, is con¬ 
sistent with our data. Note that next-nearest-neighbor 
interactions along the c axis could not be detected in the 
high-field experiment of Ref. [3TJ as their scattering plane 
was orthogonal to c. Analyzing the measured dispersion 
in DTNX, we find Jc 2 to be ferromagnetie, and therefore 
nonfrustrating. Both Jd and Jc 2 corrections are relatively 
small, and don’t qualitatively affect the general picture 
for DTN. 
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DTN 



6% DTNX 



Zapf Tsirulin 

Present 

study 

Present study Present 

study 

Present study 

et al. (2006;P3 et al. (2013^^ 

Eq. n 

only 

Eq. fl+Jd Eq. (R 

1 + Jc2 

Eq. Ill+Jd + Jc 2 



0.943 

- 

0.921 

0.923 

0.926 

0.924 

D 

0.700(3) meV 

0.767 meV 

0.682(3) meV 

0.702(2) meV 

0.702(2) meV 

0.701(3) meV 

Jc 

0.150(3) meV 

0.177 meV 

0.168(1) meV 

0.171(1) meV 

0.161(1) meV 

0.162(1) meV 

Ja 

0.015(1) meV 

0.0134(3) meV 

0.0172(5) meV 

0.0174(2) meV 

0.0167(2) meV 

0.0170(5) meV 

Jd 

- 

0.0069(9) meV 

- 

0.0064(5) meV 

- 

0.0065(4) meV 

Jc2 

- 

- 

- 

- 

-0.0103(7) meV 

-0.0104(7) meV 


D/Jo 4.7(1) 

4.34 

4.06(3) 

4.10(3) 

4.36(3) 

4.33(3) 

Ja/Jc 0.098(7) 

0.076(2) 

0.102(2) 

0.102(2) 

0.104(2) 

0.105(2) 


TABLE IL Results of fitting Eqs. ( B1][B3 ) to our 6% DTNX dispersion data and to those for the parent compound (Ref. [2^ 
Eor the parent compound, parameters deduced from the spin-wave dispersion measured in the saturated phase, are also listed.^ 


DTN 



6% DTNX 



Sizanov and 

Present 

study 

Present study Present 

study 

Present study 

Syromyatnikov (201 ij^ 

Eq. H 

only 

Eq. Q+Jd Eq. ^ 

1 + Jc2 

Eq. (|rj)+ Jd + Jc 2 


D 

0.666 meV 

0.661(1) meV 

0.661(1) meV 

0.672(2) meV 

0.672(2) meV 

Jc 

0.160 meV 

0.176(1) meV 

0.177(1) meV 

0.169(1) meV 

170(2) meV 

Ja 

0.0172 meV 

0.0190(3) meV 

0.0192(2) meV 

0.0190(3) meV 

0.0193(2) meV 

Jd 

0.0086 meV 

- 

0.0075(6) meV 

- 

0.0075(5) meV 

Jc2 

- 

- 

- 

-0.007(1) meV 

-0.007(1) meV 


D/Jo 4.16 

3.75(3) 

3.73(3) 

3.96(4) 

3.95(6) 

Ja/Jc 0.108 

0.108(2) 

0.109(2) 

0.112(2) 

0.113(3) 


TABLE III. Results of fitting the dispersion obtained through the 1/D expansion ( |B4| ) to the data for 6% DTNX. The results 
of an analogous treatment for the parent compound are also shown, following Ref. 1581 


Appendix B: Self-consistent corrections to the RPA 

The RPA result ^ can be modified to account for the 
quantum renormalization of the dispersion relatio n. In 
the framework of a generalized spin-wave approachp^^ 
relying on artificial restriction of Hilbert space for the 
spin-wave operators, Eq. (§ is transformed into: 

nuj{Q) = yD'2+4D's27(Q). (Bl) 

Structurally this formula is identical to the RPA re¬ 
sult <§• The difference is that the single-ion anisotropy 
D is replaced by the effective value D' . In addition, in¬ 
stead of the bare exchange couplings Jr, one uses renor¬ 
malized s^Jr. The values of D' and are related to 
the “true” values D and Jr by a set of self-consistent 
equations: 



BZ 


The integration is performed over the full Brillouin 
zone. Applying this correction to the parameters in Ta¬ 
ble |I] results in the values, summarized in Table |n| The 
parameters of pure DTN, determined from the spin-wave 
dispersion in the fully polarized phase,!^ are also given 
for comparison. 

An alternative approach was developed by Sizanov and 
Syromyatnikov.l^ Their dispersion equation, based upon 
the bosonization and subsequent expansion in 1/D, is: 




1 


iT)!,3p 

(27r)3 

7(P)7(K)7(P-iq,3p,3K 


BZ 


{2n) 


6 


■7(Q)-A7(Q)2 + ^^(Q)3 


li 


|7(P)^ + h(P)7(P - Q) ,3 


D2 

5 

4^ 


(27r)3 




7(P)7(K)7(P-K + Q) . . 


(27r)6 


-d-^Pd-^K. 

(B4) 

Applying this approach to our data results in the inter¬ 
action parameters, summarized in Table III Values for 
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pure DTN, as obtained in Ref. [58] by analyzing the data 
of Ref. [271 are also given for comparison. Although this 


approach leads to slightly different parameter values, the 
conclusion still holds: DTNX approaches QCP primarily 
due to the reduction of the ratio Dj 
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